Magneto-elastic waves in crystals of magnetic molecules 



O 

o 

(N 



C. Calero, E. M. Chudnovsky, and D. A. Garanin 
Department of Physics and Astronomy, 
Lehman College, City University of New York, 
250 Bedford Park Boulevard West, Bronx, New York 10468-1589, U.S.A. 

February 1, 2008 

We study magneto-elastic effects in crystals of magnetic molecules. Coupled equations of motion 
for spins and sound are derived and the possibility of strong resonant magneto-acoustic coupling 
is demonstrated. Dispersion laws for interacting linear sound and spin excitations are obtained 
for bulk and surface acoustic waves. We show that ultrasound can generate inverse population of 
spin levels. Alternatively, the decay of the inverse population of spin levels can generate ultrasound. 
Possibility of solitary waves of the magnetization accompanied by the elastic twists is demonstrated. 

PACS numbers: 75.50.Xx, 73.50.Rb, 75.45.-fj 
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I. INTRODUCTION 

Crystals of molecular magnets are paramagneta^ that 
have the ability to maintain macroscopic magnetization 
for a long time in the absence of the external magnetic 
field. This is a consequence of the magnetic bi-stability 
of individual molecules^ that in many respects behave as 
superparamagnetic particles. The latter is due to a large 
spin (e.g., 5 = 10 for Mn-12 and Fe-8) and high magnetic 
anisotropy of the molecules. Together with quantization 
of spin energy levels this leads to a distinctive feature 
of molecular magnets: A staircase hysteresis curve in a 
macroscopic measurement of the magnetizationi^ 

Hybridization of electron paramagnetic resonance 
(EPR) with longitudinal ultrasonic waves has been stud- 
ied by Jacobsen and Stevens^ within a phenomenologi- 
cal model of magneto-elastic interaction proportional to 
the magnetic field. General theory of magneto-elastic 
effects on the phonon dispersion and the sound veloc- 
ity in conventional paramagnets has been developed by 
Dohm and Fuldci^ The advantage of molecular magnets 
is that they, unlike conventional paramagnets, can be 
prepared in a variety of magnetic states even in the ab- 
sence of the magnetic field. Spontaneous transitions be- 
tween spin levels in molecular magnets are normally due 
to the emission and absorption of phonons. Interactions 
of molecular spins with phonons have been studied in the 
context of magnetic relaxation,-I i^'^'^° conservation of an- 
gular momentum, ^^'^^'^^ phonon Raman processeSj-'^ and 
phonon superradianceiJ^^ Parametric excitation of acous- 
tic modes in molecular magnets has been studied i^^d^ 
It has been suggested that surface acoustic waves can 
produce Rabi oscillations of magnetization in crystals of 
molecular magnets In this paper we study coupled dy- 
namics of paramagnetic spins and clastic deformations at 
a macroscopic level. 

When considering magneto-elastic waves in paramag- 
nets the natural question is why the adjacent spins should 
rotate in unison rather than behave independently. In 
ferromagncts the local alignment of spins is due to the 
strong exchange interaction. Due to this interaction the 
length of the local magnetization is a constant through- 



out the ferromagnet. We shall argue now that a some- 
what similar quantum effect exists in a system of weakly 
interacting two- level entities described by a fictitious spin 
1/2. Indeed, since any product of Pauli matrices reduces 
to a single Pauli matrix a a , interaction of N independent 
two-state systems with an arbitrary field A(r) should be 
linear on cTq,, 
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N 
n=l 
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where cr^") describes a two-state system located at a 
point r = r„. If A was independent of coordinates, then 
the Hamiltonian ([T|) would reduce to 



where 



N 



1/3 



r(») 
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is the total fictitious spin of the system. In this case 
the interaction Hamiltonian would commute with S^, 
thus preserving the length of the total fictitious "mag- 
netization" . This observation is crucial for understand- 
ing Dicke superradiance:^^ A system of independent two- 
state entities behaves collectively in a field whose wave- 
length significantly exceeds the size of the system. When 
the wavelength of the field is small compared to the size 
of the system but large compared to the distance be- 
tween the two-state entities, the same argument can be 
made about the rigidity of S = ^ cr^") summed up over 
the distances that are small compared to the wavelength. 
Consequently, the system that has been initially prepared 
in a state with all spins up, and then is allowed to evolve 
through interaction with a long-wave Bose field, should 
conserve the length of the local "magnetization" in the 
same way as ferromagncts do. 

The relevance of the above argument to the dynam- 
ics of magnetic molecules interacting with elastic defor- 
mations becomes obvious when only two spin levels are 
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important. This is the case when the low-energy dynam- 
ics of the molecular magnet is dominated by, e.g., tunnel 
split spin-levels or when the magneto-acoustic wave is 
generated by a pulse of sound of resonant frequency. Re- 
cently, experiments with surface acoustic waves in the 
GHz range have been performed in crystals of molecu- 
lar magnets.^'i The existing techniques, in principle, al- 
low generation of acoustic frequencies up to 100 GHz^ 
This opens the possibility of resonant interaction of gen- 
erated ultrasound with spin excitations. In this paper we 
study coupled magneto-elastic waves in the ground state 
of a crystal of molecular magnets. We derive equations 
describing macroscopic dynamics of sound and magne- 
tization and show that high-frequency ultrasound inter- 
acts strongly with molecular spins when the frequency of 
the sound equals the distance between spin levels. We 
obtain the dispersion relation for magneto-elastic waves 
and show that non-linear equations of motion also pos- 
sess solutions describing solitary waves of magnetization 
coupled to the elastic twists. 

The paper is organized as follows. The model of spin- 
phonon coupling is discussed in Section II where coupled 
magneto-elastic equation are derived. Linear magneto- 
elastic waves are studied in Section III where we obtain 
dispersion laws for bulk and surface acoustic waves. Non- 
linear solitary waves are studied in Section IV. Sugges- 
tions for experiments are made in Section V. 



II. MODEL OF MAGNETO-ELASTIC 
COUPLING 

We consider a molecular magnet interacting with 
a local crystal field described by a phenomenological 
anisotropy Hamiltonian TLa- The spin cluster is assumed 
to be more rigid than its elastic environment, so that the 
long-wave crystal deformations can only rotate it as a 
whole but cannot change its inner structure responsible 
for the parameters of the Hamiltonian "Ha ■ This approx- 
imation should apply to many molecular magnets as they 
typically have a compact magnetic core inside a large unit 
cell of the crystal. In the presence of deformations of the 
crystal lattice, given by the displacement field u(r), local 
anisotropy axes defined by the crystal field are rotated 
by the angle 

,5</,(r,i) = ^Vxu(r,i). (4) 

As a consequence of the full rotational invariance of the 
system (spins -I- crystal lattice) , the rotation of the lattice 
is equivalent to the rotation of the operator S in the 
opposite direction, which can be performed by the (25 + 
1) X (25* 4- 1) matrix in the spin space,— 

S^R-^SR, 7? = e'^*'^. (5) 

Therefore, the total Hamiltonian of a molecular magnet 
in the magnetic field B must be written as 



n 



HAe 



iS-. 



Hz + n 



ph : 



(6) 



where Ha is the anisotropy Hamiltonian in the absence 
of phonons, Ti.z — ~g^B^ ■ S is the Zceman Hamiltonian 
and Tiph is the Hamiltonian of harmonic phonons. The 
angle of rotation produced by the deformation of the lat- 
tice is small, so one can expand Hamiltonian ([6]) to first 
order in the angle 54> and obtain 



(7) 



where Ti.Q is the Hamiltonian of non-interacting spins and 
phonons 

T^o = T^s + Tiph , Ti-s = T^a + Ti-z , (8) 
and Hs^ph is the spin-phonon interaction term, given by 



-ph 



Ha, S 



(50. 



(9) 



A. Coupling of spins to the elastic twists 



For certainty, we consider a crystal of molecular mag- 
nets with the anisotropy Hamiltonian 



Ha = -DSt 



V , 



(10) 



where y is a small term that does not commute with 
the Sz operator. This term is responsible for the tunnel 
splitting. A, of the levels on resonance. 

At low temperature and small magnetic field, 
ksT^ gfisB < A, when the frequency of the displace- 
ment field u(r) satisfies u <C 2DS/h, only the two lowest 
states of Ha are involved in the evolution of the sys- 
tem. Thus, one can reduce the spin-Hamiltonian of the 
molecular magnet to an effective two-state Hamiltonian 
in terms of pseudospin-1/2 operators ct^, 



H 



1 



(11) 



where ct,; are the Pauli matrices in the basis of the Sz- 
states close to the resonance between \S) and | — S"), and 
W — Es — E^s is the energy difference for the resonant 
states at A = 0. The non-degenerate eigenfunctions of 
H^^"'^ are 



\^^) = ^{C±\S)tC^\-S)) 



with 



C+ = Wl± 



w 



(12) 



(13) 



VA2+W2 

In terms of |'0:p) the Hamiltonian (fTT|) can be written as 
-U^-lf) = _i + A2 dz , (14) 
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where ct^ are now the PauU matrices in the new basis 
|'0±), i.e., = - IV'-)(V'-|. The projection 

of the spin-phonon interaction Hamiltonian ([9]) onto this 
new two-state basis results in 

(15) 

with ay ~ + The total Hamilto- 

nian ([S]) of a single molecular magnet becomes 



7^(e//) 



-b 



H 



ph ■ 



y ■ 



(16) 



Here we have assumed that the perturbation introduced 
by the spin-phonon interaction is much smaller than the 
perturbation V producing the splitting A, which will usu- 
ally be the case. Note also that A and W can in gen- 
eral be made r-dependent to account for possible inho- 
mogeneities of the crystal. 

When considering magneto-elastic excitations we will 
need to know whether they are accompanied by a non- 
zero local magnetization of the crystal. For that rea- 
son it is important to have the magnetic moment of the 
molecule, 



(17) 



(with g being the gyromagnetic ratio and being the 
Bohr magneton), in terms of its wave function 



(18) 



where K± are arbitrary complex numbers satisfying 
\K^\'^ + \K+\'^ = 1. With the help of Eq. 1^ one obtains 



W 



A 



VW^ + A2 
A(^,) - W{a 
VW^ + A2 



{KlK- + K+K*_) 



(19) 



B. Magneto-elastic equations 

We want to describe our system of N spins in terms of 
the spin field 



AT 



n(r) = ^<5-i,5(r-r,), 



(20) 



satisiying commutation relations 

[flair), hp{r')] = 2ieQ^^n^(r)(5(r - r') . (21) 
In terms of this field the total Hamiltonian becomes 



n^-^Jd^r n(r) • b^"^^'^) (r) + K 



ph ■ (22) 



The classical pseudo-spin field n(r, t) can be defined as 

n(r,t) = (n(r)), (23) 

where (...) contains the average over quantum spin states 
and the statistical average over spins inside a small vol- 
ume around the point r. If the size of that volume is small 
compared to the wavelength of the phonon displacement 
field, then, as has been discussed in the Introduction, 
n2(r) should be approximately constant in time. Accord- 
ing to equations (fTTj) . (fTO]) and ([20|) . the magnetization 
is given by 

^ A Tij; (r) — (r) 
MJr) = 01138 ] ' — . (24) 

The dynamical equation for the classical pseudo-spin 
field n(r, t) is 



which, with the help of Eq. ([2T]l . can be written as 
.9n(r,t) 



dt 



n(r,t) X h^^ff\r,t) 



(25) 



(26) 



In this treatment we are making a common assumption 
that averaging over spin and phonon states can be done 
independently. This approximation is expected to be 
good in the long- wave limit. 

The dynamical equation for the displacement field is 



at2 



(27) 



where ctq^ = dh/deap is the stress tensor, = 
dua/dxfs is the strain tensor, h is the Hamiltonian den- 
sity of the system in 7t! = J d^rh{r), and p is the mass 
density. Note that the stress tensor has an antisymmet- 
ric part originating from the magneto-elastic interaction 
in the Hamiltonian, 



+ "a/3 ' 

is'A ny{r)ezai3 ■ 



(28) 



This implies that at each point r there is a torque per 
unit volume. 



Ta{r) = -dazSAnyir) , 



(29) 



created by the interaction with the magnetic system. 
This effect can be viewed as the local Einstein - de Haas 
effect: Spin rotation produces a torque in the crystal lat- 
tice due to the necessity to conserve angular momentum. 
With the help of equations HID, and using 

standard results of the theory of elasticity, one obtains 



V ' 

(30) 
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where ci and Ct are velocities of longitudinal and trans- 
verse sound. The source of deformation in the right hand 
side of this equation is due to the above-mentioned torque 
generated by the spin rotation. 

Equations ([26| and ([30|) describe coupled motion of 
the pseudospin field n(r,t) and the displacement field 
u(r, t). It is easy to see from these equations that in ac- 
cordance with the argument presented in the Introduc- 



tion n1 



is independent of time. It may, nev- 



ertheless, depend on coordinates, reflecting the structure 
of the initial state. In this paper we study cases in which 
the crystal of molecular magnets was initially prepared 
in the ground state n = hqCz with uq being the concen- 
tration of magnetic molecules. In this case the dynamics 
of n(r) described by equations ([26]) and (|30p reduces to 
its rotation, with the length of n(r) being a constant hq. 
Remarkably, this situation is similar to a ferromagnet, 
despite the absence of the exchange interaction. 




FIG. 1: Interacting sound and spin modes. Notice the gap 
below spin resonance ujq. 
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III. LINEAR MAGNETO-ELASTIC WAVES 
A. Bulk waves 

For magnetic molecules whose magnetic cores are more 
rigid than their environments, only the transverse part 
of the displacement field (with V • u(r) = 0) interacts 
with the magnetic degrees of freedom. This is a conse- 
quence of the fact that the elastic deformation produced 
by the rotation of local magnetization is a local twist of 
the crystal lattice, required by the conservation of angu- 
lar momentum. Let us consider then a transverse plane 
wave propagating along the X-axis. From Eqs. ([26]) and 
((50)) one obtains 

d'^u,, n d'^u,, SA dn,, 



drix 
df 
dny_ 

dt 
driz 
dt 



2p dx 



A2 



ox 



duo, 



SA 



OX 



(31) 



We shall study linear waves around the ground state 
14'+) corresponding to = nQ,nx,y = 0,Uy = 0. The 
perturbation around this state results in nonzero Ux.y 
and Uy. Linearized equations of motion are 



SA dr 



dx^ 



2p dx 



_ dn„ 



SAno^ 
ox 



h—^=-n,^/W^ + A^ . 
ot 



(32) 



(cj^ - c'^q'^)uy - iq^Uy = 
Zp 

. noSAVW^ + A2 
«9 — Uu 



Uy ~ . 



(33) 



The spectrum of coupled excitations is given by 



/ 'Z 2 2\ I i 



ft2 



In the vicinity of the resonance, 



one can write 



vW2TA2 
Ct90 = r = t^O , 



(34) 
(35) 

(36) 



with b to be determined by the dispersion relation. Sub- 
stituting equations l[35|) and p6|) into Eq. ([Mjl. one ob- 
tains 



b = ±^ 



SpCfhojo ' 



(37) 



that describes the splitting of two coupled modes at the 
resonance. The repulsion of elastic and spin modes is 
illustrated in Fig. [1] The relative splitting of the modes 
reaches maximum at W = {hujo = A): 



I npS^A _ ^ I A 



2Mc2 ' 



(38) 



where M = p/ng is the mass of the volume containing 
one molecule of spin S. Notice also another consequence 
of Eq. ([M)) : The presence of the energy gap below ujq — 



For Uy.rix.y cx exp(iga; — iwi), the above equations be- vM^^TjlA^/S (see Fig. [T]). The value of the gap follows 
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from Eq. p4|) at large q. It equals 2(5^ wq- This effect 
is qualitatively similar to the one obtained in Ref. 
from an ad hoc model of spin-phonon interaction. In 
contrast with that model our results for the splitting of 
the modes and for the gap do not contain any unknown 
interaction constants as they are uniquely determined by 
the conservation of the total angular momentum (spin + 
crystal lattice). 

According to equations p3|) and (p4)l the Fourier trans- 
forms of Uy and Uy are related through 



Uq 



UJn 



A 



quy . 



(39) 



Due to the condition of the elastic theory qUy <C 1, the 
absolute value of the ratio Uy/no is generally small, un- 
less ll! is close to luq. This means that away from the 
resonance the sound cannot significantly change the pop- 
ulation of excited spin states. At the magneto-clastic 
resonance, substituting equations (|36p and p7p into the 
above equation, one obtains: 



no 



-I \uy\ 



(40) 



Although this relation is valid only at \ny \ <^ uq, it allows 
one to estimate the amplitude of ultrasound that will 
significantly affect populations of spin states. We shall 
postpone the discussion of this effect until Section |Vl 
Meantime let us compute the magnetization generated 
by the linear elastic wave, Uy = uocos[qo{x — ctt)]), in 
resonance with our two-state spin system. The last of 
Eqs. ((32)) yields Ux — i{ijJ /toojUy. Then, with the help of 
Eq. (HH) and Eq. (gOl) one obtains 



A /2Mc? 



1/2 



qQUaCOs[q(i{x - Ctt)]. (41) 



So far we have investigated coupled magneto-elastic 
waves in the vicinity of the ground state, Uz = n^. Eqs. 
(I^Tjl also allow one to obtain the increment, F, of the 
decay of the unstable macroscopic state of the crystal, 
riz = —7^-0 J in which all molecules are initially in the ex- 
cited state IV"-)- In fact, the result can be immediately 
obtained from equations ([5^ - ([M]) by replacing no with 
—no. It is then easy to see from Eq. that in the vicin- 
ity of the resonance the frequency acquires an imaginary 
part that attains maximum at the resonance where 



to = loq{1 ± i\5\) . 



(42) 



The mode growing at the rate F = loq\5\ represents the 
decay of IV'-) spin states into spin states, separated 
by energy hw^- This decay is accompanied by the ex- 
ponential growth of the amplitude of ultrasound of fre- 
quency Wo. 



B. Surface waves 

Magneto-elastic coupling in crystals of molecular mag- 
nets can be studied with the help of surface acoustic 
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q^c 'e, 




y ; e, 



FIG. 2: Geometry of the problem with surface acoustic waves. 

waves (see Discussion). To describe the surface waves 
we chose a geometry in which the surface of interest 
is the XZ-plane and the solid extends to y > with 
waves running along the direction that makes an angle 
with the AT-axis, see Fig. [2] As usualS2, we assume 
that the displacement field u(r, t) and the components 
n2,(r, i), n^ (r, i) have the form 

A = ^oe-"fe*(«-^+«^^)e-^'^* . (43) 

It is convenient to express the components of the dis- 
placement field in the coordinate system defined by 
(e;,et,ep), see Fig.H 



Uy 



ui cos 9 — ut sin 9 

Up 

ui sin 9 + Ut cos 9 . 



(44) 



Equations of motion for u;, Ut, and Up follow from Eq. 



[uj^ + Ctia^ - q'^)]ut + — aAsin6'nj^ = 

[ 2 22 221 ■ / 2 2\ 

Lu + c^a — Ciq \ui — iaq{ci — c^ )Up 



r 2 2 2 2 21 

[w + Cia" — c^q \ Up 



— a A cos 9n,, — 
2p y 

iaq{cf - c^)ui 



Aqcos9ny — . 



(45) 



It is easy to see that for 6* 7^ fevr , k = 0, 1, 2... and Uy ^ 0, 
the transverse component ut cannot be zero, contrary 
to the case of Rayleigh waves. This is the signature of 
magneto-elastic coupling. 

As in the analysis of bulk waves, we shall study the 
linear waves around the ground state corresponding to 
the pseudospin field polarized in the Z-direction, n^ — 
no,nx,y — 0. The excitations above this state are de- 
scribed by Eqs. (|26p . which become 



— ihidUx = SA [— a(u/ cos — M( sin 0) — i(7|[ cos 0Up] 

+ i/VF2 + ny 
-ihuny = -VW'^ + A'^ Ux . (46) 

Substitution of these two equations into Eqs. (|45p leads 
to a homogeneous system of algebraic equations for w;. 
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Ut, and Up, that have a non-zero solution only if its deter- 
minant equals zero. From this condition we obtain three 
values of the coefficient a that describe the decay of the 
wave away from the surface: 



ai = 



as 




where 



(48) 



Note that if there are no spins {S = 0), then — ol2 
and one obtains decay coefficients for ordinary Rayleigh 
waves. 

The general plane wave solution for the components of 
the displacement field can be written as 



Ey(fc)g-Qfcygj(<3^2;+(j^z)g-iwt ^ 

fe=l 



(49) 



i*i = Ui(x — vt), Ui = riiix — vt). For such a choice, Eq. 
(ED) gives 



du. 
dx 



V 



SA 



2p(c2 - z;2) 



(51) 



where x = x — vt and the constant of integration was put 
zero assuming that there is no duy/dx independent from 
Substituting this into the equations of motion for n, 



(47) Eqs. dni), one obtains 



dUy 

driz 



where 



;VTy2 ^ A2 
hv 



52A2 



The system of Eqs. ((5^ can be reduced to 



(52) 



2p(c2 - t-2)Viy2_^^2 



(53) 



where is the amplitude corresponding to each ak 
and i — l,t,p. For each k, the amplitudes wJq ■* , u|g\ Upg"* 
are related through Eqs. ([l5|) (there are two independent 
equations, so we can express, e.g., u^Q^Upg'' in terms of 

(k) \ 

). Therefore, there still are three unknowns, say 

„(2) ,,(3) 



''10 ' "iO ' ^10 



The boundary conditions for the stress 
tensor at the surface, aiy\y=o = 0, provide a system of 

homogeneous equations for u'iq' , u|g' and u|q' , whose de- 
terminant must be zero to allow for non-trivial solution. 
From this last condition we obtain the dispersion relation 
for surface magneto-elastic waves: 




OJ 



3/2 



/ 9252^2^0 C0S2 e + 2Mc2;i(w2 _ ^2)(g2 _ ^2/^2) 



2,2 



q'^ sin^ ( 



0. 



(50) 



This equation should be solved numerically to obtain the 
dispersion law for magneto-elastic modes. Qualitatively, 
the repulsion of the modes is similar to the one shown in 

Fig.m 



IV. NON-LINEAR MAGNETO-ELASTIC 
WAVES 

An interesting feature of Eqs. (pij) is the existence of 
transverse non-linear plane wave solutions of the form 



d Uy 



I 2 2 

7C -7 Hy 



(54) 
(55) 



where C is a constant of integration. The first integral 
of the last differential equation is 



dn., 

di 



^[^) =-:,{l-lC)nl-:Ln^y + A>0, (56) 



where A is another integration constant. 

We are interested in real bounded solutions of Eq. ((55|) 
with Uy vanishing at x—vt — > ±00, so that the integration 
constant A must be zero. In this case, for the right hand 
side of Eq. ([55)1 to be positive we must have 1 — 7C < 0. 
Then, the solution of Eq. ((55|) is 



From the equations 
dn 



(57) 



(58) 



one determines with the help of the condition n'^ + riy + 
n\ = Hq that C = ±no . Therefore, 7 must satisfy 
I7I > 1 /no for the equation ([55]) to have a solution satisfy- 
ing the conditions specified above. Setting the reference 
point ^0 = one obtains 



nyiO = ±|^Vl7K-l sech [a/ItK-I C 



(59) 



7 



M^/ (gAiBSno) 



V. DISCUSSION 




FIG. 3: Magnetization inside the soliton as a function of ^ for 
W = 0. 



SO that 



n,(0 = ±lT2 



\l\ 



— sech^ 



- 1 d • (60) 



In these formulas, the upper sign corresponds to 7 > 
and the lower sign to 7 < 0. 

Eq. ((60)) describes a solitary wave of a characteristic 
width 



^0 



I 



hv 



(61) 



travelling at a speed v. The parameter 7 given by Eq. 
()53|) is determined by f , which is the only free parameter 
of the soliton. The magnetization inside the soliton is 
given by Eq. (j24|) with Ux and defined by equations 
dSHl) - dSQl). At, e.g., = 

25(|7|no-l) 



xsech 



tanh 



The condition 



A 



> 1 



(62) 



(63) 



requires v to be very close to the speed of sound Ct ■ This 
is a consequence of A being very small compared to Mcf . 
Note that the maximal value of the magnetization inside 
the soliton. 



\Mz\ = gfJ-sS no 



l7l 



(64) 



is, in general, of the order of saturation magnetization 
Mq = g^sSnQ. We should also note that although the 
above non-linear solution of the equations of motion for- 
mally allows V to be both slightly lower or slightly higher 
than Cf, the supersonic soliton should be unstable with 
respect to Cherenkov radiation of sound waves. 



Eq. (|38)) provides the splitting of the bulk sound fre- 
quency in a magnetized crystal of magnetic molecules in 
the vicinity of the resonance between sound and spin lev- 
els. At a zero field bias {W = 0) the resonant condition, 
A = Hctq, should be easily accessible at low A. However, 
the splitting given by Eq. ([55)) will be very small unless A 
is in the GHz range or higher. Such a large A will be also 
beneficial for decreasing inhomogeneous broadening of A 
and for insuring low decoherence of quantum spin states. 
Surface acoustic waves can, in principle, be generated up 
to lOOGHz^^. They may also be easier to use for the 
observation of the discussed splitting. By order of mag- 
nitude it will still be given by Eq. ([38)) . Substituting into 
this equation S = 10, A ~ O.I K (frequency / in the GHz 
range), Mc^ ~ lO^K, one obtains Smax 10^^. This will 
be observable if the quality factor of ultrasound in the 
GHz range exceeds 100. The magneto-elastic nature of 
the splitting can be confirmed through its dependence on 
the angle between the wave vector and the easy magne- 
tization axis of the crystal, see Sec. III-B. Observation 
of the gap, 2d'^ujo, in the excitation spectrum (see Fig. 
(T)) will be more challenging. For practical values of S the 
gap is likely to be small compared to the width of the 
spin resonance and the width of the ultrasonic mode in 
the GHz range. 

Eq. (gO)) shows that at M - lO^^ig ^j^^^ _ lO^^s'^ 
ultrasound of amplitude uq ~ 0.1 nm will significantly 
affect population of spin levels. Moreover, it will result 
in the oscillating magnetization of large amplitude, Eq. 
([¥!)) . We have also demonstrated that one can prepare 
the crystal in the excited spin state and generate ultra- 
sound due to the decay of the population of that state. 
This result is another confirmation of the phonon laser 
effect suggested in Ref. ilSj. Equations ([^ and ([55)) show 
that at uiQ ^ lO^^s""'^ the amplitude of the sound wave 
may grow at the rate as high as F ~ lO^s^^. Magneto- 
elastic effects studied in this paper should be sensitive 
to the decoherence of spin states. However, when the 
oscillation of spin population is driven by the external 
acoustic wave, the latter should force the phase coherence 
upon the spin system. To provide the resonance condi- 
tion, the broadening of the level splitting due to disorder 
and dipolar fields should be small compared to A. If it 
is not, the tunnel splitting, A, should be increased by 
applying a sufficiently large transverse magnetic field. 

One fascinating prediction of our theory is the exis- 
tence in molecular magnets of solitary waves of the mag- 
netization reversal coupled to elastic twists. Such waves 
have quantum origin as they are related to the quantum 
splitting of spin-up and spin-down states. They can be 
ignited in experiment that starts with all molecules in the 
ground state. Such a state of the crystal has zero magne- 
tization as the molecules are in a superposition of spin- 
up and spin-down states. The soliton discussed above 
is characterized by a narrow region of a large non-zero 
magnetization that propagates through the solid with the 
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velocity close to the speed of transverse sound. It can be 
generated by, e.g., a localized pulse of the magnetic field 
or by a localized mechanical twist, and detected through 
local measurements of the magnetization. In general the 
width of the soliton, given by Eq. ((6T|) . is of order of the 
wavelength of sound of frequency yW^ + A^/ h, though 
wider solitons are allowed if |7|rto — *■ 1. In experiment 
this width should depend on the width of the field pulse 
or the size of the twisted region that generates the soli- 



ton. 
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